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ABSTRACT 

Local and intermediate redshift {z ~ 0.5) galaxy samples obey well correlated 
relations between the stellar population luminosity and maximal galaxy rotation that 
define the Tully-Fisher (TF) relation. Consensus is starting to be reached on the TF 
relation at z ~ 0.5, but work at significantly higher redshifts is even more challeng- 
ing, and has been limited by small galaxy sample sizes, the intrinsic scatter of galaxy 
properties, and increasing observational uncertainties. We present here the TF mea- 
surements of 41 galaxies at relatively high redshift, spectroscopically observed with the 
Keck/DEIMOS instrument by the DEEP2 project, a survey which will eventually offer 
a large galaxy sample of the greatest depth and number yet achieved towards this pur- 
pose. The 'first-look' sample analyzed here has a redshift range of 0.75 < z < 1.3 with 
(z) — 0.85 and an intrinsic magnitude range from Mb of —22.66 to —20.57 (Vega). 
We find that compared to local fiducial samples, a brightening of 1.5 magnitudes is 
observed, and consistent with passive evolutionary models. 

Key words: galaxies: evolution - galaxies: formation - galaxies: kinematics and 
dynamics - galaxies: spiral 



1 INTRODUCTION 



The obse rved evolution of the Tully-Fisher relation (TF; 
ITuUv fc~F isher 1977) has the potential to provide consider- 
able insight into how the properties of today's galaxy popu- 
lation were established. In particular, this relation between 
the luminosity and rotation velocity of disk galaxies - probed 
through the total mass-to-light ratio of galaxies - reveals 
links between dark matter haloes and the stellar populations 
which form within them. This should ideally allow us to dis- 
tinguish between the various proposed theories of galaxy 
formation and evolution. 

In the hierarchical structure formation paradigm of cold 
dark matter (CDM), small galactic haloes (10^"^" Mq) first 
condensed out at peaks in the near-uniform dark matter dis- 
tribution. Eventually, through sufficient numbers of mergers, 
these sub-units accumulated to produce the massive haloes 
inhabited by the ~ L* galaxies that dominate the stellar 
mass density from 2; ~ 1 to the present. One possibility 
is that these first, s mall haloes were the sites of significant 
early star formation ||Mo et al.lll998l : ISommer-Larsen et al.l 
|2003| ). and the resulting stellar mass subsequently accumu- 
lated through mergers (along with their dark matter haloes) 
into more massive galaxies without the need for substantial 



later star-formation. In this case, the growth of stellar mass 
in individual galaxies was driven primarily by mergers of 
already luminous sub-units, and hence the observed mass- 
to-light ratio should have been relatively steady over time, 
evolving only slowly due to the opposing actions of passive 
aging of the stellar populations and residual star-formation. 

In contrast, there is some growing evidence that star 
for mation was sup pressed in smaller haloes at early times 
( Nelan et al. 1 l2005|). The suggested cause s of this range 
from the increasing ultraviolet background (IWvithe fc Loeb 
200£ ), the effects of star formation feedback (jSpringel et al. 
20051 ) (particularly from the theoretically expected initial 
population of very-massive, zero-metallicity stars), or even 
AGN prefer entially en hancing star formation in more mas- 
sive haloes l|Silkll2005D . Whichever mechanism is responsi- 
ble, this implies that galaxies must have formed a substan- 
tial amount of their stellar mass afterwards, at later times, 
following the many mergers that assembled most of their 
dark matter haloes. This is somewhat reminiscent of the 
old monolithic coll apse model for massive galaxy formation 
l|Eggen et al.lll963 ). and implies substantial evolution in the 
mass-to-light ratio of galaxies - as stars form in the already 
assembled massive halo. This would also require more star 



2 K. Chiu, S.P. Bamford, A. Bunker 



formation per unit mass in smaller haloes than in more mas- 
sive ones at later times, the so called 'downsizing' effect. 

The relative order of mass assembly and star formation 
is cleary of significant interest. Therefore, examining evolu- 
tion in the mass-to-light ratio of galaxies has been pursued as 
an important indicator of the galaxy formation mechanism. 
Such evolution is often probed through observed changes in 
fundamental scaling laws versus redshift, such as the Faber- 
Jackson and Kormendy relations for el lipticals (believed to 
be pr ojections of a fundamental plane, IDiorgovski fc Davis! 
Il987l'). For spiral gal axies, the TuUy-Fisher relation (TF) 
i Tullv fc Fished Il977l ) is the primary scaling law. Statisti- 
cally well-populated TF diagrams over a range of redshifts 
should be able to test the various galaxy formation options 
through changes in the slope, intercept, and scatter of the 
fitted TF relation. 

TuUy-Fisher relations have been relatively well char- 
acterized for g alaxies nearby, originally for distance indi- 
cator purposes (iPierce fc TullvlllQQ^ : iGiovanelli et al.lll997l : 
iDale et al.ll200ll '). but recently in a more inclusive manner for 
assessing galaxy evol ution, and to compare with higher red- 
shift work (e.g.. iKannappan et al.. 2002, '). A number of stud- 
ies have measured the TF relation at intermediate redshifts 
(z ~ 0.4-1) and attempted to infer its evolution with cos- 
mic time. However, a consensus on the TF relation evolution 
over the past Hubble time (out to z ~ 1) has not yet been 
achiev ed, although very recent work such as IWeiner et al.l 
l|2006! ) have greatly improved upon the few existing high 
redshift studies. Difficulty arises due to the relatively low 
numbers of galaxies typically used; the intrinsic scatter in 
the TF relation means that many tens of galaxies are needed 
per redshift bin for a reliable measurement, particularly to 
constrain the slope of the relation. The primary cause of dis- 
crepancies between studies appears to be differing selection 
criteria, and the manner in which these have been corrected 
for, if at all. The applied corrections can vary widely depend- 
ing on the percieved aim of each particular study. In addi- 
tion, the use of different local relations, TF fitting methods, 
and internal extinction corrections has complicated compar- 
isons between studies. 

As a result, various groups have identified significant 
evolution i n the i3-ba nd TF relation o ut to z 1 (e.g., 
Vogt et all [19 96 ■ 1997: iRix et al.! 1 19971: [Simard fc Pritchet 



Harden ct al. 2003; B ohm et al.ll2004l : lB6hm fc Ziegler 
200&: .Nakamura et al. 20oH )^ or state an upper limit 



l|Bamford et al.! 12006! ) on the evolution, while others find 
no evolution once selectio n effects have been accounted for 
llSimard fc Pritchet! [1999I : IVogt 2001 ) or kine matically dis- 
turbed galaxies discarded (jPlores et al.![200^ ). However, a 
consistent picture is begining to emerge, as these studies 
become more mature. One fairly certain aspect is that there 
has been rather little evolution in the intercept of the TF, 
ranging from no evolution to a modest brightening of ~ 1 
mag at fixed rotation velocity by z ~ 1. This is much less 
than would be expected if the star formation rates of mas- 
sive disk galaxies have evolved as strongly as is measured 
for t he global star formation rate density of the universe 
ifBamford ct al. 2006 ). 

Constraining evolution in the TF slope is more dif- 
ficult, though again any evolution is fairly modest. Any 
overall brightening of the TF that is observed appears 
to be mainly driven by low mass galaxies, as noticed by 
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Figure 1. Redshift distribution of galaxies measured in this work 
(solid line). Dotted line indicates distribution of full DEEP2 spec- 
troscopic sample, prior to TF analysis selection, discussed in text. 
The full sample distribution has been scaled by 0.05 for compar- 
ison purposes. 



iBamford" et al.l ||2006! ). identified as an additional popula- 
tion by IVogtl d 200lh. and r eflecte d in the TF slope chang e 
measured by iBohm et all (|2004l '): iBohm fc Zieglei! (|2006! ). 
This is particularly the case when considered along with 
the observed lack of evol ution in the i^-band a nd stellar 
mass TF relations found bv lConselice et al.1 (|2005| ). implying 
that while the stellar to total mass of galaxies stays roughly 
constant, the colour, and hence luminosity- weighted age or 
fraction of stars that have recently formed, has evolved, pos- 
sibly with a dependence on galaxy mass. However, it is still 
unclear how much of this effect is actually due to differen- 
tial luminosity evolution at a given rotation velocity, or a 
symptom of underestimating the rotation velocities for an 
increasing population of kinematically disturbed galaxies at 
hig her redshifts. 

iFlores et all (|2006! l use integral field spectroscopy to in- 
vestigate this question, and also find that the Tf-band, or 
stellar mass, TF relation has not evolved for a small sam- 
ple of undisturbed disks out to z ~ 0.7, while galaxies with 
signs of kinematical disturbances are generally measured at 
low rotation velocities for their K-hanA luminosity. In the 
B-band, however, they show some evidence for a brighten- 
ing of 0.5 mag, a nd a tentative chan ge of slope consistent 
with the studies o f iBohm et al. I (|2004! ). even for undisturbed 
disks. [Smith et ahl (2004) find a consistent result from near- 
infrared integral field spectroscopy for a galaxy at z = 0.82. 

A novel method in TF studies is the use of strong 
lensing by galaxy clu sters to magnify distant galaxies 
(|Swinbank et all 12003! ). This enables the examination of 
z ~ 1 galaxy rotati on curves with r esolut ion comparable 
to z ~ 0.1 studies. ISwinbaiik et al.! (I2006D use this tech- 
nique to study six z ~ 1 galaxies, four of which are found to 
have regular rotation curves. These are found to define an 
7-band TF relation very similar to that found locally, and 
a B-band TF with similar slope to the local relation, but 
offset to brighter magnitudes by ~ 0.4-0.5 mag. 
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Figure 2. Image and spectral model fitting results for the galaxies in the present TF sample are shown, with DEEP2 galaxy identification 
names. The first 3 columns illustrate a) the original galaxy postage stamp image, b) the fitted model using an exponential disk profile in 
GALFIT, and c) the model-subtracted image. The last two columns show d) the original spectral flux at the 3727A [Oil] line, followed 
by e) the model generated by the ELFIT2PY procedure. Galaxy images are 60 X 60 pixels (12" X 12"), while the 2d spectra are 80 pixels 
wide (9", or 25A). 



While such integral field work is important in providing 
a complete map of the galaxy velocity field (an improve- 
ment on the spatial losses accompanying the usual tech- 
nique of long-slit spectroscopy) the number of galaxies ob- 
servable with this method is rather limited by present in- 
strumental capabilities. To measure rotation curves for a 
large sample of galaxies requires spatially-resolved multi- 
object spectroscopy, and currently the most efficient method 
is the use of slit-masks with individual slits tilted to align 
with the galaxy m ajor axes. The DE IMQS optical spectro- 
graph on Keck II (|Faber et al.ll2003l ) provides this capabil- 
ity, and a major pro ject exploiting thi s instrument has been 
the DEEP2 project l|Davis et al.ll2004l '). This survey has un- 
dertaken spatially-resolved multi-slit spectroscopy of a large 
sample of galaxies colour selected to lie at 0.7 < z < 1.4. 
The red sensitivity of the instrument, coupled with its sta- 
bility and lack of fringing, means that [Oil] 3727 A rotation 
curves can be traced cleanly throughout this redshift range. 

As well as spectroscopy, in order to construct a TF rela- 



tion at intermediate redshift one requires high quality imag- 
ing, to provide magnitudes and importantly for measuring 
the inclination angles of the targetted galaxies, needed to 
correct the rotation velocities for projection effects. 

In this paper, we measure the rotation curves of 41 
galaxies at (z) ~ 0.85 from the DRl release of the DEEP2 
spectroscopy, and use them to determine the evolution of 
the B-band TF relation. In Section 2 we discuss the DEEP2 
spectroscopy from Keck/DEIMOS, and the imaging from 
CFHT. In Section 3 we detail the fitting of galaxies in the 
imaging in order to derive inclinations, and the emission- 
line modelling to measure the rotation velocities. Section 4 
presents our resulting TuUy-Fisher relation and a compari- 
son with previous work at low and high redshift. Our con- 
clusions are in Section 5. In this paper all magnitudes are 
expressed based on Vega normalization, and throughout we 
adopt the standard 'concordance' cosmology of f^A/ = 0.3, 
f^A = 0.7, Ho = 70 km s"^ Mpc"\ 
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Figure 3. - continued. 



2 DATA SOURCES AND PROCESSING 

In this work, we aim to compare the fitted TF parame- 
ters for our present sample with those for similar samples 
at intermediate and lower redshifts. The main uncertain- 
ties in constructing a TF relation usually arise from the 
galaxy rotation velocities. These are determined from the 
rotation curves, typically using the shapes of emission lines 
measured in the spatially-resolved two-dimensional spec- 
troscopy. These velocities are then corrected for projection 
effects due to the inclination of the galactic disk, and mis- 
alignment of the slit from the major axis of the galaxy, both 
determined from morphological fits to the imaging data. 
Therefore, while the intrinsic scatter of the galaxies them- 
selves place limits on the uncertainties of the fitted TF rela- 
tion, galaxy imaging and spectroscopic data of high signal- 
to-noise and resolution should mi nimize the remaining fun - 
damental sources of uncertainty (|Kannappan et al]|2002l ). 
In particular, in this work we take advantage of the high 
spatial (0'.'l2/pixel) and wavelength resolution (68kms~^) 



of the DEEP2 survey, combined with the wide area imag- 
ing archive over the DEEP2 survey region from the CFHT 
CFH12K mosaic imager in the /-band (corresponding to the 
rest-frame B-band at z ~ 0.9). 



2.1 DEEP2 spectroscopy 

We began our sample selection by examining the DEEP2 
dataset for suitable galaxies for TF measurement. The 
DEEP2 spectroscopic observations were based on an optical 
imaging survey taken from 1999-200 using the CFH12K 
mosaic camera l|CuiUandre et al.ll200ll ) in the BRI bands to 
a depth of Rab < 24.1. Galaxy target selection and slitmask 
design were conducted using this imaging (which covered 
four fields composing a total of 3 deg^), with color cuts in 
{B ~ R) & {R~ I) to optimi ze the selection o f galxies in the 
redshift range z = 0.75-1.4 (|Coil et al.ll2"ooi ). 

Although the DEEP2 project will eventually obtain 
DEIMOS spectroscopy for some 40,000 galaxies over the 
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total area described above, the current data release con- 
tains ~ 7500 galaxies over ^ 1.2 deg^. Of these, approx- 
imately 4300 have confident redshifts as measured by the 
DEEP2 pipeline software (Q = 4), and we inspected these 
galaxies for strong emission lines at the expected position 
of [OII]3727A in each spectrum. In this first ("IHS") phase 
of the DEEP2 project, one hour spectroscopic integrations 
were obtained on each 16' x 4' slitmask, containing about 
85 objects per mask. Later phases of DEEP2 ("3HS") will 
integrate on a smaller galaxy sample for > 3 hours. After 
being dispersed by the 1200/mm grating, the wavelength 
coverage is 6500 - 9100 A, so [OII]3727A is detectable over 
0.77 < z < 1.42, nearly the entire redshift range targeted 
by the BRI colour-cuts. The spectra have a typical spectral 
resolution of i? = A/AAfwhm ~ 4000 ~ 68kms~^ deter- 
mined by the 1" slit size, sampled by ~ 0.3A pixels. The 
spatial pixel size is 0'.'119. Of the 4300 spectra, 388 were 
found to have strong [Oil] emission in the aforementioned 
wavelength range, and suitable for measuring rotation ve- 
locities using the line fitting task with integrated single- line 
signal of S/N > 20. 

The spectra were provided by the DEEP2 team in the 
DRl public releas^il as a fully reduced dataset, with catalogs 
and wavelength-calibrated multi-extension FITS 2d-spectra 
and Id spectral extractions. We did little additional pro- 
cessing except that prior to the model fitting procedures de- 
scribed below, we took the additional step of rectifying the 
2d spectra to have orthogonally aligned spatial and spectral 
axes, with constant pixel scale. This was carried out to cor- 
rect spatial distortions and because the slits of the DEIMOS 
instrument can be tilted to accomodate the position angles 
of the selected galaxies (within ~ 30°, constrained by the 
CCD columns and desired spatial resolution). This results 
in a generally diagonal (or even curved) wavelength solution 
across the spatial rows. Using the DEEP2 wavelength so- 
lution at the given spatial center of each galaxy in the slit 
mask, the solution was reproduced in all rows, and the fiux 
interpolated and rebinned to a constant dispersion of 0.32 
A/pixel. 

2.2 CFH12K imaging 

With the sample of suitable spectra identified, we extracted 
images of the galaxies for measurement of the necessary 
shape parameters. Searching the CFHT archive for the three 
currently released DEEP2 spectroscopic fields (out of four 
total to be made public), we identified useful imaging expo- 
sures taken by the CFH12K instrument. The CFH12K cam- 
era has a platescale of O'.'206/pixel, and covers a 42' x 28' 
rectangular field of view using a 2 x 6 array of 4Kx2K 
CCDs, for a total of 0.326 square degrees per pointing. Two 
fields, "0230-1-0000" and "2330+0000" , had /-band imaging 
which overlapped with most of the DEEP2 spectra where 
we identified strong [Oil] emission. With a median redshift 
of (z) = 0.85 for the galaxy sample used here (Figure [T|, 
the choice of the /-band allows the photometry to be eas- 
ily converted to the i3-band rest-frame absolute magnitude 
commonly used for TF studies, with minimal /("-correction 
uncertainties. 

^ http://deep.berkeley.edu/DRl 



The individual 6 x 10 minute dithered imaging expo- 
sures in each field were combined using standard image re- 
duction techniques, including the use of bad pixel masks 
and calibration frames ("master detrend" bias/dark/flat- 
field frames) identified for each camera run in the CFH12K 
archive. Cosmic ray rejection was implemented in addition 
to bad pixel masking and local sky subtraction during the 
median combining of the separate exposures. Finally, each 
reduced imag e was astrometrica lly corrected to the USNO- 
Bl.O catalog (|Monet et aLlbOOSl ). using the WCSTools rou- 
tines of D. Mink, which automatically update the appropri- 
ate header WCS information given a standard star catalog 
and chosen reference stars in the image of interest. This 
allowed automated extraction of the image of each galaxy 
selected this work, based on the J2000 RA/Dec positions 
given in the DEEP2 spectroscopic catalogs. At this point, 
290 (out of 388) galaxies with strong [Oil] fiux were also cov- 
ered by these two CFH12K imaging mosaics, and remained 
in the sample for further analysis. 



3 DATA ANALYSIS 

For the 290 objects selected to have strong [Oil] line emis- 
sion in the DEIMOS spectroscopy, and with CFH12K /- 
band images, we used the imaging to model the morpho- 
logical parameters (giving the disk inclincation angle, i, and 
position angle) , and the spectroscopy to model the emission- 
line rotation curve (yielding Kot sini). We describe both of 
these below. 

The TF relation also requires absolute magnitudes {Mb 
in this instance), and we determined these from the /-band 
imaging, introducing a redshift-dependent /iT-correction to 
match the restframe-ZJ with the observed-/ filter. Finally 
we corrected the absolute magnitude for intrinsic extinction 
in the disk galaxy, depending on inclination and luminosity. 

3.1 Galaxy image and spectral profile fitting 
procedures 

First, image shape modeling wa s carried ou t on ea ch galaxy 
using the GALFIT program of iPeng et al] l|2002l ). in order 
to calculate the position angle (PA) and axial ratio values 
required in subsequent steps. For each galaxy, a 100 x 100 
pixel cutout (12" x 12") was made from the CFH12K image, 
preserving astrometric and flux information. A point spread 
function (PSF) was produced by stacking ~ 10 good point 
sources within each CFII12K array chip, resulting in a profile 
with a FWHM of 2.5 pixels (0.52"). 

GALFIT convolves this PSF with a model galaxy im- 
age based on input initial parameter estimates {x, y posi- 
tion, effective radius, magnitude, axial ratio, position an- 
gle) and perturbs these parameters until the residuals of 
the model-subtracted image are minimized. The initial in- 

! iut parameters were det ermined by the SExtractor program 
Bertin fc Arnoutslll996l ) on the galaxy images. An exponen- 
tial disk profile was chosen as the desired galaxy model, and 
the program iterated to completion, returning each of the 
minimized input parameters with errors. Of the 290 galax- 
ies input to GALFIT, half returned acceptable fits, while 
the remaining galaxies generally failed to converge due to 
unsalvagable fiux profiles (i.e., galaxies not well described 



6 K. Chiu, S.P. Bamford, A. Bunker 



by exponential disks, or unresolved in the image). Exam- 
ples of galaxy profile fits, and the raw/subtracted CFH12K 
images are shown in Figure [2] 

The fitting of faint emission lines in order to extract re- 
liable rotation velocities is a challenging task, and highlights 
the desirability of high resolution and high signal-to-noise 2- 
dimensional spectra in order to determine confident galaxy 
rotation curves. In order to extract as much information as 
possible from the emission lines, and to consistently model 
the effects of inclination, seeing, finite slit width and in- 
strumental resolution, we employed an enhariced version of 
the E LFIT2PY routine developed by SPB (jBamford ct al. 
120051 ). which is b ased partly on the synthet ic emission line 
fitting software of lSimard fc Pritchetl ()l998h . 

The software is given the galaxy inclination, seeing, in- 
strumental spectral resolution, slit position angle misalign- 
ment with respect to the galaxy major axis, slit tilt in 
the mask, and produces synthetic emission lines by con- 
structing model 2D galaxy flux and velocity profiles, as- 
suming a set of parameters, and then mimicking the ob- 
servations. ELFIT2PY uses the Metropolis searching algo- 
rithm to find the set of parameters which minimize the 
residuals between the synthetic and observed emission lines. 
The searched parameters are the rotation velocity, the emis- 
sion flux scalelength (assuming an exponential disk), the 
turnover radius of the rotation curve, and the spatial and 
wavelength centres of the emission line. The use of more 
complicated flux proflles was attempted, but it was found 
that these generally could not be constrained. An intrin- 
sic rotation curv e of the form Vrot/{r°' + r^Y^"" is assumed 
l|Courteaulll"997l '). with a = 5 (found to best reproduce the 
observed rotation c urve shapes, thoug h ill-constrained, and 
in agreement with iBohm et aLll2004l l. and where r is the 
radius and rt is the turnover radius. Such an intrinsic ro- 
tation curve is necessary in this case, rather than the dis- 
continuous form used previously with ELFIT2PY, in order 
to adquately fit the high resolution of the DEEP2 spectra. 
Note that the emission flux scalelength and the turnover 
radius of the intrinsic rotation curve were allowed to vary 
independently, which was frequently found to be necessary. 
The total flux in the line was measured by integrating the 
background- and continuum-subtracted image of the line, 
and held fixed for the fit, with a constant [Oil] doublet ra- 
tio of /([Oil] 3727.092 A)//([OII] 3729.875 A) = 0.8. The 
background level was also fixed at zero. 

In previous versions of ELFIT2PY, the chi-squared 
comparison between the synthetic and observed lines was 
done on a pixel by pixel basis in the sp ectra. However, some 
other groups (e.g., iBohm et al.l |2003 ) first trace the syn- 
thetic and observed emission lines, by fitting the position of 
the line in each spatial row independently, and then com- 
pare these traces to produce a chi-squared value used to 
judge the goodness of fit. This has the advantage of remov- 
ing one level of dependence on the assumed galaxy flux pro- 
file, though it still of course plays a role in the construction 
of the synthetic line. It also weights the fit to preferentially 
reproduce the observed rotation curve rather than the flux 
profile, which may be affected by discrete star-formation re- 
gions and asymmetries. However, this is at the cost of leaving 
the emission scalelength unconstrained, which must there- 
fore by fixed relative to the rotation curve turnover radius. 
Therefore a hybrid approch was adopted to take advantage 



of the high resolution of the DEEP2 spectra (which allow a 
high quality trace of the emission line to be measured) , while 
somewhat reducing the influence of asymmetries in the ob- 
served flux proflle which are apparently more prevalent at 
the high redshifts of the current sample. 

This hybrid method combined the chi-squared values 
from pixel-to-pixel comparisons of the synthetic and ob- 
served spectra, and those from comparisons of both the trace 
centres and FWHM of the synthetic and observed lines. The 
variation of the observed FWHM in each spatial row helps 
constrain the rotation curve as it is a measure of the gra- 
dient of the rotation curve at each spatial position. If the 
rotation velocity is increasing quickly with radius, the ob- 
served FWHM in that spatial row is higher. This combined 
chi-squared value was the quantity used by the Metropolis 
algorithm to determine the best-fltting parameters and their 
confidence intervals. 

The traces were performed by fitting a double Gaussian 
function to each spatial row independently. The Gaussians 
had a fixed separation appropriate for the [OH] doublet at 
the redshift of the galaxy, and the same FWHM. Only the 
centre of the doublet, single line FWHM and amplitudes of 
the two components were allowed to vary. Fits were judged 
to be reliable if the Levenberg-Marquardt least-squares min- 
imisation converged from the appropriately chosen initial 
parameters, if the errors on the best-fitting centre and am- 
plitudes were sensible, and if the best-fitting FWHM was 
larger than the spectral resolution. 

Note that the velocity convention in this work is to use 
the Kot "half-velocity" term, which specifies the maximum 
absolute rotation velocity, of the asymptotic assumed intrin- 
sic rotation curve, with respect to the center of the galaxy, 
rather than the velocity- width, 2Kot, traditionally used to 
express the total velocity gradient measured across the entire 
galaxy. Examples of synthetic emission-line images, created 
using the best-fit parameters, appear in Figure [21 

The emission lines of the 290 galaxies selected thus far 
were inspected visually in more detail before fitting, and 
those with no sign of being spatially extended were re- 
jected. Of the remaining galaxies, those with unsuccessful 
GALFIT measurements of their inclination and position an- 
gle were also rejected. Then, only galaxies with inclinations 
more than 30° from face-on, and hence a reasonable com- 
ponent of the rotation velocity along the line-of-sight, were 
fit with ELFIT2PY. In addition, galaxies with a large mis- 
alignment angle (> 30°) between their major axes and the 
slit were rejected, to avoid large and uncertain corrections 
for this effect. 78 galaxies were thus remaining to be fit with 
ELFIT2PY. 

In a number of cases, the best fitting parameters from 
ELFIT2PY indicated that the rotation velocity was unre- 
liable, for example due to the emission fiux profile being 
consistent with a point source, or the end of the reliably 
measurable emission line occuring at smaller radii than the 
best-fitting intrinsic rotation curve turnover radius. In or- 
der to remove such uncertain galaxies from the sample we 
imposed the following quality criteria. 

If the scalelength of the fitted emission fiux exponential 
disk profile is too small, then the fit is unlikely to be useful, 
and thus is rejected. The limit adopted is 0.14" - the angle 
subtended by one pixel for a galaxy observed with a slit tilt 
of 30°, the worst case we admit. This corresponds to 1-1.2 
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Table 1. Extracted galaxy parameters. Tabulated for each galaxy: name, coordinates in J2000, redshift, apparent magnitude, absolute magni- 
tude, galaxy inclination (i), angular difference between galaxy PA and spectroscopic slit orientation, angular difference between slit mask and 

spcctroscoy)ir slit, iiiaxiiiial rotation \'olocity from SA'stoiiiic. 



Galaxy ID 


RA 


Dec 


Redshift 




nil 


Mb 




i 


Slit-Gal 


Mask-Slit 


Vrot 




(J2000) 


(J2000) 






(Vega) 


(Vega) 




n 


n 


n 


(km/s) 


32001201 


23:29:40.78 


-00:01:41.9 


0.83580 


22, 


.43 + 0.03 


-20.85 + 0. 


11 


55.94 + 11.1 


7.4 


28.0 


79.8 + 7.93 


32001370 


23:29:23.38 


-00:01:33.8 


1.15568 


21 


.52 + 0.01 


-22.66 + 0. 


15 


50.20 + 8.90 


13.4 


-16.4 


44.2 + 3.31 


32001373 


23:29:20.31 


-00:00:58.4 


1.03207 


22 


.15 ± 0.01 


-21.74 + 0. 


37 


55.24 + 31.6 


12.5 


-30.0 


225.3 + 54.07 


32001527 


23:29:13.43 


-00:01:04.3 


1.11611 


21 


.93 + 0.01 


-22.40 + 0, 


14 


61.31 + 6.50 


9.6 


6.8 


156.3 + 5.94 


32001538 


23:29:05.50 


-00:00:55.0 


1.05244 


22, 


.72 ± 0.02 


-20.94 ± 0. 


,17 


40.53 + 25.8 


14.8 


30.0 


148.3 + 51.18 


32001550 


23:29:17.52 


-00:01:42.9 


1.05390 


22 


.35 ± 0.03 


-21.99 + 0. 


,13 


68.89 ± 4.90 


14.9 


0.0 


164.0 ± 3.14 


32001607 


23:28:51.26 


-00:01:08.9 


1.05507 


22 


.54 ± 0.02 


-21.38 + 0. 


,13 


55.24 ± 9.70 


4.2 


9.8 


175.1 ± 10.86 


32015354 


23:30:02.10 


+00:08:17.3 


0.83198 


21 


.98 ± 0.01 


-21.11 ±0. 


09 


43.94 ± 13.3 


0.0 


-30.0 


198.5 + 25.75 


32015753 


23:29:50.20 


+00:06:47.1 


0.82388 


21 


.71 +0.01 


-21.51 + 0, 


11 


52.41 + 11.6 


-0.7 


-30.0 


361.0 + 29.29 


32015948 


23:29:55.86 


+00:07:25.7 


0.95589 


22 


.35 ± 0.02 


-21.50 + 0, 


,28 


63.89 + 17.9 


-16.1 


-30.0 


138.2 + 11.45 


32016140 


23:29:36.71 


+00:08:56.1 


0.80110 


21 


.91 ±0.01 


-21.57 + 0, 


10 


65.79 + 6.20 


-25.4 


-30.0 


85.8 + 6.53 


32016292 


23:29:40.37 


+00:07:21.5 


0.83904 


22 


.49 ± 0.02 


-21.05 + 0, 


15 


66.42 + 10.0 


4.0 


26.7 


124.9 + 11.75 


32016833 


23:29:30.78 


+00:07:56.4 


0.78703 


22 


.25 + 0.03 


-21.47 + 0, 


18 


74.33 + 7.10 


6.1 


-30.0 


252.8 + 34.23 


32017515 


23:28:57.81 


+00:08:41.9 


0.99842 


22, 


.02 ± 0.01 


-21.76 + 0. 


,11 


55.24 + 8.30 


-21.6 


30.0 


187.3 + 9.72 


32020384 


23:30:19.50 


+00:11:06.8 


1.24909 


22 


.63 ± 0.02 


-22.26 ± 0. 


,36 


68.89 ± 14.7 


29.3 


-13.2 


170.8 ± 12.97 


32021318 


23:29:56.02 


+00:09:41.6 


0.80140 


22 


.57 + 0.02 


-20.57 ± 0. 


,16 


44.76 ± 25.2 


20.8 


0.0 


65 ± 19.00 


42001530 


02:30:43.59 


+00:22:49.2 


0.75225 


21 


.89 ± 0.01 


-21.12 + 0. 


,52 


52.41 + 56.1 


24.0 


0.0 


166.9 + 158.20 


42001854 


02:30:32.79 


+00:21:53.5 


0.80186 


21 


.51 + 0.01 


-21.49 + 0, 


,25 


42.26 + 37.3 


7.0 


18.5 


157.8 + 99.05 


42001899 


02:30:24.87 


+00:21:36.3 


0.90477 


21 


.93 + 0.01 


-22.15 + 0, 


,66 


72.54 + 23.0 


18.2 


30.0 


163.6 + 11.36 


42002654 


02:30:02.52 


+00:22:49.7 


0.80614 


21 


.46 + 0.01 


-21.91 + 0, 


13 


60.00 + 9.20 


2.6 


20.7 


215.6 + 10.64 


42006919 


02:30:40.80 


+00:24:27.8 


0.89.347 


21 


.23 + 0.01 


-22.33 + 0, 


10 


55.94 + 8.30 


7.1 


-30.0 


182.4 + 9.12 


42007501 


02:30:30.40 


+00:24:45.1 


0.92616 


21 


.67 + 0.01 


-22.01 + 0, 


12 


58.66 + 9.40 


-4.4 


18.3 


109.1 + 6.92 


42013571 


02:30:50.69 


+00:29:17.0 


0.82416 


22, 


.10 + 0.01 


-20.96 + 0. 


,07 


43.11 + 10.1 


12.2 


3.4 


214.6 + 20.95 


42014180 


02:30:34.44 


+00:28:07.2 


0.75354 


21 


.74 ± 0.01 


-21.22 + 0. 


,17 


49.45 ± 19.9 


-25.6 


0.0 


161.3 + 26.94 


42014642 


02:30:27.62 


+00:29:19.9 


0.75022 


22 


.50 ± 0.02 


-20.74 ± 0. 


,62 


65.79 + 36.1 


16.1 


0.0 


150.3 ± 26.90 


42015874 


02:30:02.02 


+00:29:18.1 


0.86724 


22 


.52 + 0.02 


-20.81 ± 0. 


,50 


54.54 + 51.0 


-14.6 


17.4 


218.9 ± 133.10 


42016553 


02:29:48.27 


+00:29:39.6 


0.80765 


21 


.59 + 0.01 


-21.56 + 0. 


,25 


50.20 ± 27.6 


-18.9 


0.0 


139.4 ± 34.85 


42016993 


02:29:32.82 


+00:29:31.9 


0.88262 


22, 


.41 ± 0.02 


-21.17 + 0. 


31 


62.61 +22.1 


-23.0 


-30.0 


64.2 + 19.36 


42017441 


02:29:17.55 


+00:28:27.9 


0.85701 


22 


.49 + 0.02 


-21.,32 + 0, 


,46 


72.54 + 19.3 


-0.8 


4.2 


138.0 + 14.36 


42017561 


02:29:13.87 


+00:27:12.9 


0.85383 


22 


.26 + 0.02 


-21.01 + 0, 


,40 


52.41 + 43.6 


-6.1 


0.0 


69.4 + 33.44 


42023413 


02:29:59.75 


+00:33:20.1 


0.81719 


21 


.46 + 0.01 


-21.89 + 0, 


,24 


57.99 + 19.1 


20.1 


30.0 


179.9 + 20.20 


42024446 


02:29:28.61 


+00:32:01.6 


0.81243 


21, 


.04 ± 0.00 


-22.06 + 0. 


,09 


44.76 + 11.4 


-3.6 


30.0 


210.1+22.09 


42027940 


02:31:01.50 


+00:35:10.5 


0.86831 


22 


.39 ± 0.02 


-21.46 ± 0. 


,17 


45.57 ± 24.8 


-2.4 


30.0 


216.6 + 55.86 


42035205 


02:31:09.36 


+00:40:21.8 


0.86450 


22 


.76 + 0.02 


-20.58 ± 0. 


,17 


55.94 + 18.1 


14.4 


7.7 


134.7 ± 15.69 


42036185 


02:30:59.17 


+00:40:03.4 


0.79288 


22 


.53 + 0.02 


-20.82 + 0. 


,48 


59.33 ± 40.4 


25.3 


-30.0 


203.1 +59.06 


42036908 


02:30:45.18 


+00:38:43.5 


1.29168 


22 


.90 ± 0.02 


-21.84 + 0. 


,32 


62.61 ± 16.8 


1.1 


0.0 


164.9 ± 13.84 


42040167 


02:29:28.90 


+00:37:14.8 


0.78000 


21 


.26 ± 0.01 


-22.09 + 0. 


,16 


61.31 + 10.4 


26.5 


30.0 


171.1 + 8.80 


42043879 


02:30:48.83 


+00:43:05.4 


0.77335 


21 


.46 ± 0.00 


-21.71 + 0. 


15 


55.94 + 12.5 


-6.0 


11.1 


184.1 + 14.67 


42044594 


02:30:37.44 


+00:42:10.3 


0.90493 


21 


.50 + 0.00 


-22.48 + 0. 


15 


68.89 + 7.30 


-1.5 


-12.2 


181.8 + 4.58 


42053518 


02:30:07.45 


+00:44:46.2 


0.92744 


21 


.72 + 0.01 


-21.74 + 0. 


,22 


48.70 + 25.0 


9.8 


30.0 


102.8 + 23.77 


42058430 


02:30:57.37 


+00:48:03.1 


0.97522 


21, 


.83 + 0.01 


-22.05 + 0. 


,52 


60.65 + 33.8 


17.0 


-30.0 


186.5 + 38.66 



kpc for the redshifts of our objects. In order to also reject 
galaxies with uncertain emission scalelengths which are close 
this limit, we reject any galaxy for which the fit scalelength 
is consistent with being below 0.14" at the Icr level. 

In order to consider how far out in each galaxy we mea- 
sure the emission line to, we define the extent of the line 
as half of the spatial distance over which the line is reliably 
traced, after removing single 'good' points with no 'good' 
neighbours as spurious. If the extent is less than that one 
would expect for a point source given the seeing, which wc 
take as 0.5" then we reject that galaxy from further consid- 
eration. 

Finally, we also consider whether the transition from 
the inner rotation curve slope to the flat region is observed 
with sufficient signal-to-noise. This is done by requiring that 
the radius at which the fit intrinsic rotation curve turns over 



is within the reliable extent of the line. Furthermore, in the 
case where the turnover is not observed, the turnover radius 
will be highly uncertain. To exclude such galaxies we require 
that the turnover radius is inconsistent with being beyond 
the lino extent with at least la confidence. 

Following the application of these quality criteria, 41 
galaxies remained with reliable rotation velocity measure- 
ments. 

3.2 Galaxy photometry, if-correction, and 
extinction correction 

In addition to the rotation velocity, the second input re- 
quired for the TF relation is the galaxy absolute magni- 
tude. In order to compare with other samples, while avoid- 
ing large if -corrections to the photometry, we use the ab- 
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solute magnitude in the rest-frame _B-band, Mb, in the AB 
zero-point system. For our galaxies, with median a redshift 
of {z) — 0.85, this roughly corresponds to the observed- 
frame I filter (centered at 8100 A). For each galaxy of 
interest, the 7-band photometry was extracted from the 
DEEP2 photo metric catalogs accompanying the spectro- 
scopic dataset (|Coil et al.ll2004 ). These catalogs were pro- 
duced using th e IMCAT im aging and photometric reduc- 
tion software of lKaiser et al.i (1995 ) , having been finally cal- 
ibrated to the SDSS standard system. As an independent 
check on the 7-band photometry, we also verified magni- 
tudes during our SExtractor and GALFIT proceedures. 

Using sy nthetic template spectra of Sbc, Scd, and Im 
galaxy types (jBruzual fc Charlotll200 j ). appropriate for the 
range of disk galaxies with star formation and hence [Oil] 
emission, we simulated flux through the rest-frame B{z = 0) 
and redshifted I{z) filters. We applied the average resulting 
7f-correction and error to each galaxy according to its red- 
shift. 

Finally, we applied an inclination-dependent extinction 
correction to ea c h gal axy. Historically, t he co rrections of 
iTuUv fc Fouau3 (Il985l ) and iTullv et all (|l998l ) have been 
used most commonly, the latter incorporating an additional 
dependen ce on galaxy lumin osity (or rotation velocity) . We 
adopt the lTullv et"an (|l998l ) method, for comparison of our 
work to recent TF samples at similar a nd lo wer redshifts. 
The extinction is defined in lTuUv etal] (|l998l ) to be: 



with 



A = 'yB log(a/6). 



7s = -0.35(15.3H- Ms) 



where a/b is the galaxy major-to-minor axial ratio. Table [T] 
displays all of the measured parameters of the galaxies in 
our sample. 



4 THE TULLY-FISHER RELATION 

From the original catalog of 4300 galaxies with high quality 
redshifts in the DEEP2 DRl data release, the aforemen- 
tioned steps have reduced the sample significantly, specif- 
ically through inspection of the 2d spectra for [Oil] emis- 
sion (N — 400), imaging data overlap matching (A'^ — 290), 
galaxy image profile modeling {N = 150), and suitable in- 
clinations and slit misalignments for emission line modeling 
(A'^ = 78). Following rejection of those galaxies for which 
we judge the measured rotation velocity to be unreliable, 
as described in section 13.11 we now have a useful sample 
of 41 galaxies with full photometry, morphology and ro- 
tation information. The sample spans a redshift range of 
0.75 < z < 1.3 with (2) — 0.85 and an intrinsic Mb magni- 
tude range of -22.66 to -20.57. 

From this sample, we carried out several least-squares 
fittings to a straight line function: Mb = a + blogVrot, 
where a and b are the intercept and slope of the TF re- 
lation, respectively. The slope and/or intercept value, and 
scatt er width were minim ized as the free parameters, and 
as m iBamford et"all (120061 ) . we use weights of 1 /a^ for each 
point, where = crf^gy^^^ -I- b^ffM^ + crfnt, and ai„t is the 
intrinsic scatter of the TF relation. One data point was re- 
jected by sigma-clipping, specifically object 32001370, due 




log(V^„^ [km/s]) 

Figure 4. The TF relation of the 41 galaxies in the present 
sample are plotted in absolute Ms magnitude versus maximal 
rotation velocity, and overlaid with the resulting linear TF fits 
discussed in text (in order of decreasing steepness, solid red line: 
free fit, dotted blue line: slope of VOl, dashed green line: slope 
of TPOO). Light dotted blue and dashed green lines indicate the 
original local VOl and TPOO relations, again in order of decreasing 
steepness. Object 32001370 (extreme upper left) was excluded 
from fits as a > 3(T outlier. 



to its unusual position more than 3a from the main sam- 
ple. Three fits were then carried out - first, a free-fit of all 
parameters was allowed. Because this yielded a relatively 
poorly constrained slope, two subsequent fits were made 
with slop e fixed to local reference samples, specifically those 
of TuUy fc Pierce! (|2000l l or lVerhei icn ( 2001). E ach of these 
compa rison samples employed a version of the ITuIIv et al.l 
(| 19981 ) extinction correction, as is used for our sample here. 
Table 2 lists these fitted parameters, and Figure 3] displays 
the resulting TF distribution and fits. 

These fitted parameters present several interesting as- 
pects in comparison to the local fiducial TF relation. First, 
an overall brightening of 1.5 ±0.2 mag is seen in this sample 
at (z) ~ 0.8 compared to z '- -j galaxies at pre sent. This 
brightening is consistent with iBohm et al. I (|2004|) who find 
an off set of —(1.22 ± 0.56 mag)2;, as well as Bamford et al.l 
(|2006t ). who find an offset of (-1.0±0.5 mag)2. The strength 
of the effect may be explained in part by a colour selection 
effect in the DEEP2 galaxy sample, which leads to observa- 
tion of relatively younger stellar population galaxies (and in 
turn, decreased mass-to- light ratios). A larger sample from 
the full DEEP2 d ata may allow this colour- TF relation to 
be evaluated, as in^Kannapp an fc BartonI (|2004l ). Also of in- 
terest is that the brightening of this sample compared to the 
local fiducial samples seems to be strongest at the bright end 
of the TF r elation (Mb < -20 .5). This is in contrast with 
the result of lBohm et al.l (|2004h . in which little brightening 
(~ 0.5 mag) is found the relatively bright galaxies, but sig- 
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Table 2. Best-fit TF values for the sample of 41 DEEP2 galaxies. Parameters below fit the relation Mb = a + b(log Vrot +c), discussed 
in text, where c is an additional offset factor to account for the use of log Vi-ot versus log 21^0* in TPOO or VOl. 





a 


b 


c 


total scatter 


intrinsic scatter 










mags (dex) 


mags (dex) 


Free-fit 


-21.57 ±11.67 


-12.06 ±6.51 


-2.2 


1.72 (0.143) 


1.52 (0.126) 


TPOO 


-21.58 ±0.17 


-7.27 ±0.00 


-2.2 


1.07 (0.147) 


0.91 (0.125) 


vol 


-21.57 ±0.20 


-9.00 ±0.00 


-2.2 


1.30 (0.144) 


1.12 (0.124) 



nificant brightening found in the faint end galaxies 1.5 
mag, Mb ~ —18). 

Does this magnitude offset accord with the hierarchi- 
cal picture, in which small galaxies with pre-formed stellar 
populations merge into progressively more massive systems, 
whose mass-to-light ratios change only slowly with the dim- 
ming of the passively evolving galaxies? We have calculated 
the possible contribution of such passive evolution to ob- 
served changes in high redshift galaxy luminosities and find 
that such dimming over time is consistent. Using Bruzual 
& Chariot exponential decay models, we find that young 
galaxies at 2 = 0.85 (with luminosity-weighted ages of 3-5 
Gyr) will have dimmed by AAfs ~ 1.0 mag over the 6.5 
Gyr to 2 = 0. Hence we conclude that the observed bright- 
ening of the TF relation in this sample is consistent with 
pure luminosity evolution. 

We additionally find that the amount of the brighten- 
ing offset is within the intrinsic scatter of the sample, and 
comparable to the scatter found in the local reference popu- 
lations discussed above. Because intrinsic scatter in the TF 
relation reflects variations in the mass-to-light ratio of in- 
dividual galaxies, stellar mass fraction, and deviations from 
ordered rotation, these effects may be useful to separate in 
larger similar samples. 

A few points are worth consideration in the future. 
First, it will be useful to examine whether TF samples are 
really dominated by disk galaxies, and the extent to which 
unidentified irregulars may be influencing and scattering the 
relation. This may be partially influenced by the selection 
of the present sample, which is based on the color cut pro- 
ducing the DEEP2 galaxies as well as the use of relatively 
bright galaxies in order to accomplish the spectral and spa- 
tial fltting tasks here. Also, while this and other TF studies 
examine the blue, star-formation dominated region of the 
galaxy spectrum, the redder portion of galaxy flux displays 
tighter correlation and may be more useful for detection of 
parameter evolution. The effects of momentarily bright star 
formation episodes (mostly influencing the restframe blue 
bands) could be mitigated. Such work (at 2 = 1 for exam- 
ple) would benefit from H or Ji'-band imaging, correspond- 
ing to I band in the re stframe, as has been demonstrated by 
IConselice et"aLl l|2005h . Finally, with future improvements in 
spectroscopic and imaging spatial resolution (such as with 
HST), it may be possible to more thoroughly separate the 
disk and bulge flux components of examined galaxies, and 
even address non-uniform star-forming mass regions. 



5 CONCLUSIONS 

The population of 41 reliable galaxies examined here is a sig- 
nificant addition to the existing sample at (2) ^ 0.85, with 



spectral resolution of J? ~ 4000 allowing the determination 
of the TF parameters and their evolution. Our detection of 
a 1.5 mag brightening offset in the TF intercept compared 
to the pres ent 2 = relation is consistent wit h the previ- 
ous work of lBohm et all ||2004^■ iBamford et all (2006), and 
the recent results of Weiner et al.l 1 20061 ) . Wesuggest that 
this is consistent with passive evolution of galactic sub-units 
that have already undergone star-formation prior to having 
assembled into the large galaxies we see today. This is in con- 
trast to the scenario where significant new star formation is 
proposed to occur within the dark halos that assembled first. 

Even if galaxies experience only short episodes of star 
formation activity, and then evolve passively without further 
mergers, we have demonstrated that changes in correlated 
parameters, such as in the TF relation, should be observable. 
As stellar populations age, their decreasing UV luminosities 
ought to be reflected through changes in galactic mass-to- 
light ratios; conversely, activity through mergers and rekin- 
dled stellar activity leaves similarly observable signatures. 

In the future, it may be possible to separate larger sam- 
ples of similar galaxies into redshift bins spanning the entire 
high redshift range of interest, in order to further quantify 
the evolution of the TF intercept. The DEEP2 survey pro- 
vides a promising dataset to carry out such work, and to 
extend such studies towards examining the role of individ- 
ual galaxy variations in the overall scatter of TF and similar 
relations. And while our and other concurrent studies are 
capable of revealing evolution of the TF linear fit intercept 
versus redshift, in principle, changes in the TF fit slope it- 
self could also be probed more accurately than has yet been 
possible - perhaps revealing differences in useful observables 
such as star formation efficiency versus time. 
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